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Abstract. The electronic transport properties in the presence of a temperature gradient 
VT in disordered systems near the metal-insulator transition (MIT) are considered. The 
d.c. conductivity a, the thermoelectric power S, the thermal conductivity K and the Lorenz 
number Lq are calculated for the three-dimensional (3D) Anderson model of localization 
using the Chester-Thellung-Kubo-Greenwood formulation of linear response. We show that 
a, S, K and Lo can be scaled to one-parameter scaling curves with a single scaling parameter 
kBT/\{^- E,)/E^\. 
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1 Introduction 

In this paper we shall demonstrate that in the Anderson model of localization near the 
metal-insulator transition (MIT) the d.c. conductivity cr, the thermoelectric power S, 
the thermal conductivity K and the Lorenz number Lq obey scaling. 

A fundamental model often used in dealing with transport phenomena of disordered 
media is given by the Anderson Hamiltonian jl]. Investigations of this model yield 
that the energy E dependent conductivity a{E) behaves as 



1 - 



\E\ < E, 
0, 1^1 > E, 



in energy regions near the mobility edge Ec at which the MIT occurs g]. Here ctq is 
a constant and is a universal critical exponent. Using Eq. (|l]) we have been able to 
study the temperature T dependence of cr, S, K and near the MIT |^ using linear 
response theory B. Here we present our observations of their scaling features. 



2 Calculating the transport properties 

In the presence of a small temperature gradient VT in an open circuit, S is the coef- 
ficient of proportionality between VT and the electric field it induces. The coefficient 
K relates VT to the heat current while the Lorenz number Tq measures the ratio be- 
tween K and the product uT . The transport properties are defined in the framework 
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of linear response theory as 
a — Lii, 

L22L11 — L21L12 



S 



K 



and L 







\e\TLn ' 
L22L11 — L21L12 



(2) 



where e is the electron charge and ks is Boltzmann's constant. For a noninteracting 
system with no inelastic scattering processes present, such as phonon mediated hop- 
ping, the kinetic coefficients are given in the Chester-Thellung-Kubo-Greenwood 
formulation Ml as 



Uj = (-1)*+^' / A{E) [E y^{T)] 



i+j-2 



dI{E,ii,T) 
dE 



dE , 



(3) 



where i,j = 1,2, /x is the chemical potential of the system, f{E,iJ,,T) is the Fermi 
distribution function, and A{E) describes the system dependent features. In the An- 
derson model, one finds A{E) to be equal to the critical behavior (P of a near the 
MIT II . 

Previously, we have numerically determined /i(r) for the 3D Anderson model from 
its density of states Using that result, i.e., /i(T) oc T^, it is straightforward to 
determine the kinetic coefficients and, consequently, the transport properties (|^). 

In order to compare our results with experiments we have chosen the energy units 
to be in electron volts. For the same reason, we measure a in units of fl^^cm^^ . In 
this paper we only show the results obtained for disorder W = 12 where W is the 
width of the box distribution of randomly chosen potential energies in the Anderson 
model of localization. We have obtained similar results for other disorders not too 
close to the critical disorder Wc ~ 16.5 |^ and where there are no large fluctuations 
in the density of states. For W = 12, the corresponding mobility edge is at Ec = 7.5 

The value of is set to 1.3 in agreement with numerical results Note that this 
choice of v determines the magnitude of the transport properties but does not change 
their T dependence 1^ upon which their scaling is based. 



3 Scaling Features 

Wegner |^ was the first to show for disordered noninteracting electron systems that a 
can be written close to the MIT in a scaling form as 

a{t, n) = b^-'^F{b^/''t, ¥n) . (4) 

Here -F is a system dependent function of the dimensionless parameter t giving the 
distance from the critical point, is an external parameter such as the frequency, b 
is a scaling parameter, d is the dimension and z is the dynamical exponent. For the 
noninteracting case, z ^ d [p^ . In this paper, the appropriate parameters are d = 3, 
t=\{^j.-Ec)/Ec\,fl = T and b = t-" and Eq. (|) simplifies to 
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Fig. 1 Left: The d.c. conductivity obtained in the metallic regime: Ef ~ Ec = —0.01 eV 
(o), -0.007 eV (v), -0.005 eV (□), -0.003 eV (A), -0.001 eV (O), in the critical regime: 
Ef - Ec = eV (•), and in the insulating regime: Ef - Ec = 0.001 eV (+), 0.003 eV (x). 
The thin solid line is the scaling curve onto which all data from the different energy regimes 
collapses. The thick solid line denotes a high temperature expansion ^ while the dashed 
line is the Sommerfeld expansion [|, ^ for Ef- Ec^ -0.01 eV. Right: The scaling plot of 
the thermoelectric power obtained from the different regions in the same way as for a. 

In Fig. we show that the a data from the metaUic {\Ef\ < Ec), the insulating 
{\Ep\ > Ec) and the critical {Ep = Ec) energy regions collapse onto a single scaling 
curve when plotted as a function of kBT/\{fi — Ec)/ Ec\- In the right hand side of Fig. 
^ we illustrate similar scaling of S. Note that there is no pre-factor in the scaling form 
such as in cr since S becomes independent of T at the MIT as T — > 0. It was pointed 
out in Rcf. 1^ that S for high and low T might scale in terms of ksT /\Ep — Ec\- This 
is approximately valid |^ although the proper scaling variable should be kgT/KiJ. — 
Ec)/Ec\- The scaling curve for S converges at the value 228.4 //V/K as predicted in 
Ref. for J/ = 1.3. In the metallic regime, our result agrees with the Sommerfeld 
expansion result at low T. We find that S grows infinitely large in the insulating 
regime. This can be attributed to the decrease in charge carriers which contribute to 
a current. 

The corresponding scaling curves for K and Lq are shown in Fig. ||. They also 
scale with respect to ksT/Kid — Ec)/Ec\- The normalization factor in K is due to 
the fact that K oc ksT"^^ as T — > at the critical regime. Lq does not require a 
pre-factor since it is independent of T at the MIT ||ll| . The result for Lq shows that 
it approaches 7r^/3 in the metallic regime as expected in the Sommerfeld free electron 
theory of metals. In the critical regime, Lq converges to 2.414 which is the value for 
= 1.3 predicted in Ref. jl^. Due to the exponential decay of the derivative of the 
Fermi function as T — > 0, Lq approaches v + 1 = 2.3 in the insulating regime. 

To summarize, we achieve scaling for the four transport quantities cr, S, K and 
Lq with respect to a single parameter ksT/Kfj, — Ec)/Ec\- This factor stems from 
the form of Ly and from the fact that the explicit T dependence enters via fj, and 
the Fermi distribution function. Comparing our results from a in Fig. I with Eq. 
(||) we obtain vz = 1. With v = 1.3 this yields a dynamical exponent z — 0.77 
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Fig. 2 Left: Scaling plot of the thermal conductivity obtained from the same energy regions 
as given in Fig. ^ The thin solid line is the scaling curve onto which all data from the 
different energy regimes collapses. The thick solid line denotes a high temperature expansion 
j^, Isj while the dashed line is the Sommerfeld expansion ^ for Ef — Ec = —0.01 eV. 
Right: The scaling of the Lorenz number. The additional symbol (*) denotes data obtained 
deeper in the insulating regime: Ef — Ec = 0.01 eV. 



which is to be compared to z = d = 3 as expected from scaling arguments [lO| . 
Moreover, our resuh that vz — 1 seems to violate the Harris criterion jisj which 
requires I'z > 2. However, this observation needs to be examined in more detail since 
the Harris criterion might be altered for our choice of T dependence via the Fermi 
distribution. 
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